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ABSTRACT 

We compare the observed mass functions and age distributions of star clusters 
in six well-studied galaxies: the Milky Way, Magellanic Clouds, M83, M51, and 
Antennae. In combination, these distributions span wide ranges of mass and age: 
10 2 < M/M < 10 6 and 10 6 < r/yr < 10 9 . We confirm that the distributions are 
well represented by power laws: dN/dM oc M@ with /3 ~ —1.9 and dN/dr oc r 7 
with 7 rs —0.8. The mass and age distributions are approximately independent of 
each other, ruling out simple models of mass- dependent disruption. As expected, 
there are minor differences among the exponents, at a level close to the true 
uncertainties, ep ~ e 7 ~ 0.1-0.2. However, the overwhelming impression is the 
similarity of the mass functions and age distributions of clusters in these different 
galaxies, including giant and dwarf, quiescent and interacting galaxies. This is 
an important empirical result, justifying terms such as "universal" or "quasi- 
universal." We provide a partial theoretical explanation for these observations in 
terms of physical processes operating during the formation and disruption of the 
clusters, including star formation and feedback, subsequent stellar mass loss, and 
tidal interactions with passing molecular clouds. A full explanation will require 
additional information about the molecular clumps and star clusters in galaxies 
beyond the Milky Way. 

Subject headings: galaxies: individual (Magellanic Clouds, Milky Way, M51, 
M83, NGC 4038/39) — galaxies: star clusters: general — stars: formation 



1. INTRODUCTION 

Star clusters form in the dense parts ("clumps") of molecular clouds (Lada & Lada 2003; 
McKee & Ostriker 2007). They are subsequently destroyed by several processes, beginning 
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with the expulsion of residual gas by massive young stars ("feedback"), further mass loss 
from intermediate- and low-mass stars, tidal disturbances by passing molecular clouds, and 
stellar escape driven by internal two-body relaxation (Binney & Tremaine 2008). These 
processes disperse the stars from clusters into the surrounding stellar field. Star clusters 
therefore represent an intermediate stage in the transformation of the clumpy interstellar 
medium (ISM) of a galaxy into a relatively smooth stellar distribution. 

In this paper, we show that there are some remarkable similarities in the statistical prop- 
erties of cluster populations in different galaxies. In particular, we focus on the univariate 
mass and age distributions, ip(M) oc dN/dM and x( T ) dN/dr, and the bivariate mass-age 
distribution, g(M,r) oc d 2 N/dMdr. These distributions encode valuable information about 
the formation and disruption of clusters, for example, whether less massive ones dissolve 
faster than more massive ones. In the next section, we present the mass and age distribu- 
tions of the clusters in six well-studied galaxies, based on published data, but now displayed 
in a uniform manner. In the following section, we provide a partial theoretical explanation 
for these observations, and in the last section, we summarize and place our results in the 
context of other work in this field. 

As in our previous papers, we use the term "cluster" for any concentrated aggregate of 
stars with a density much higher than that of the surrounding stellar field, whether or not 
it also contains gas, and whether or not it is gravitationally bound. This is the standard 
definition in the star formation community (Lada & Lada 2003; McKee & Ostriker 2007). 
Some authors use the term "cluster" only for gas-free or gravitationally bound objects. We 
reject these definitions for several reasons. (1) They ignore the fact that gas-free and gas- 
rich clusters (including molecular clumps and HII regions) are the same objects observed in 
different evolutionary phases. (2) It is impossible to tell from images alone which clusters 
are gravitationally bound (have negative energy) and which are unbound (have positive 
energy). In principle, spectra would help, but in practice, they are seldom available and 
are often contaminated by non-virial motions (stellar winds, binary stars, etc). (3) iV- 
body simulations show that unbound clusters retain the appearance of bound clusters for 
many (~10-50) crossing times (Baumgardt & Kroupa 2007). By "disruption," we mean the 
removal of mass from a cluster, whether this occurs gradually or suddenly, and whether it 
leaves the cluster bound or unbound. 

2. MASS AND AGE DISTRIBUTIONS 

Figures 1 and 2 show the mass functions and age distributions of the clusters in six 
nearby galaxies. The mass functions are plotted for several disjoint age intervals and the age 
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distributions for several disjoint mass intervals. These figures are based on data from Lada 
& Lada (2003) for the solar neighborhood in the Milky Way, Chandar et al. (2010a) for the 
Large and Small Magellanic Clouds (LMC and SMC), Chandar et al. (2010c) for a large 
field in M83, Chandar et al. (2011) for M51, and Fall et al. (2009) for the Antennae. For 
the clusters in the solar neighborhood, the masses and ages were derived from photometry 
of individual stars; for the clusters in the other galaxies, they were derived from integrated 
photometry in several wavebands (usually UBVIHa) and comparisons with stellar popu- 
lation models. The Ha photometry is crucial for distinguishing clusters younger and older 
than ~10 7 yr. Lada & Lada (2003) and Chandar et al. (2010b) describe these methods in 
more detail. 

For ease of comparison, we have made two simple adjustments to the published mass 
and age distributions when constructing Figures 1 and 2. First, we replotted them in a 
uniform format: log(dN/dM) against log(M/M ) and \og(dN/dr) against log(r/yr). For the 
solar neighborhood, the original distributions were presented in the form log(M dN/d log M) 
against log(M/M ) and log(c?A^/(ilogr) against log(r/yr). Second, we adopted a uniform 
conversion from light to mass based on stellar population models with the Chabrier (2003) 
IMF. For the LMC, SMC, M51, and Antennae, the original distributions were based on 
models with the Salpeter (1955) IMF, which have AlogM = 0.2 and Alogr = 0.0 relative 
to models with the Chabrier (2003) IMF. 

The observed mass and age distributions are well represented by featureless power laws: 



We list the best-fit exponents and their formal la errors for the 12 mass functions and 10 age 
distributions in Tables 1 and 2. The straight lines in Figures 1 and 2 show the corresponding 
power laws. Both the mean and median exponents for this sample are /3 = —1.9 and 
7 = —0.8, and the standard deviations of individual exponents about the means are ap = 0.15 
and a 7 = 0.18 (with full ranges -2.24 < < -1.70 and -1.05 < 7 < -0.54). As a result 
of stochastic fluctuations in the luminosities and colors of clusters, the true uncertainties 
(errors) in the exponents, ep and e 7 , are usually larger than the formal la errors listed in 
Tables 1 and 2, with typical values ~ e 7 ~ 0.1-0.2 (Fouesneau et al. 2012)0 Since these 
are similar to the dispersions ap and <r 7 , we cannot tell whether the small differences among 
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allowance for these effects, the true uncertainties increase to tp ~ e 7 ~ 0.2. 
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the exponents are real, although we do expect differences at roughly this level, as explained 
below. 

Figures 1 and 2 show that the mass and age distributions are essentially independent of 
each other. This follows from the parallelism of the mass functions in different age intervals 
and the age distributions in different mass intervals. The vertical spacing between the age 
distributions differs only because the adopted mass intervals differ, a consequence of the 
different distances, limiting magnitudes, and sample sizes among the galaxies. Thus, we 
can approximate the bivariate mass-age distribution as a product of the two univariate 
distributions: 

g{M, r) oc i/j(M) X (t) oc MV. (3) 

We originally introduced this model for clusters in the Antennae, with the speculation that 
it might also apply to clusters in other galaxies (Fall et al. 2005; Fall 2006; Whitmore et al. 
2007). We emphasize that, for each galaxy, our power-law fits are based on data in large, but 
limited, ranges of mass and age, typically AlogM ~ Alogr ~ 2. In combination, however, 
they span even wider ranges, roughly 10 2 < M/M & < 10 6 and 10 6 < r/yr < 10 9 . 

The age distribution, in general, reflects the difference between the formation and disrup- 
tion rates of clusters. We know from observations that the star formation rate has remained 
roughly constant, within factors of 2-3, during the past 10 9 yr in the solar neighborhood 
(Wyse 2009 and references therein), LMC (Harris & Zaritsky 2009), and SMC (Harris & 
Zaritsky 2004). Simulations indicate that this rate has also varied slowly in the Antennae 
for at least the past 10 8 yr and possibly 10 9 yr (Karl et al. 2011). The formation rates of 
stars and clusters likely track each other closely, simply because most stars form in clusters 
(Lada & Lada 2003; McKee & Ostriker 2007). A factor of 2-3 variation in the formation rate 
of clusters between r ~ 10 6 yr and 10 9 yr would change the exponent of the age distribution 
by only A7 0.10-0.17. This is similar to the dispersion among galaxies and the true 
uncertainty mentioned above (<r 7 ~ e 7 ~ 0.1-0.2) and is much smaller than the magnitude 
of the typical exponent, I7I « 0.8. Thus, we are confident that most of the observed decline 
in the age distributions is caused by the disruption of clusters rather than variations in their 
formation rate. 

Two important conclusions now follow directly from Figures 1 and 2 and Equations (1)- 
(3). First, the shape of the mass function ijj(M) is preserved over wide ranges of age, although 
its amplitude, proportional to the age distribution x( t )j declines by a large factor. Thus, 
there is no evidence for mass-dependent disruption of clusters in this sample of well-studied 
galaxies. Second, both the mass and age distributions are remarkably similar from one galaxy 
to another, with differences mostly attributable to observational uncertainties. Thus, there 
is little if any evidence for galaxy-dependent disruption of clusters in this sample, which 
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includes giant and dwarf, quiescent and interacting galaxies. This degree of similarity in an 
astronomical setting probably justifies the adjective "universal." However, since we expect 
and possibly observe minor differences among the mass and age distributions in different 
galaxies, we sometimes employ the more modest term "quasi-universal." 

It is also worth asking what the mass and age distributions would look like if the 
disruption were mass dependent. Figure 3 helps to answer this question. The curves in this 
diagram are based on a model in which clusters form at a constant rate with a power-law 
initial mass function and are then disrupted gradually at the following mass-dependent rate: 



Boutloukos & Lamers (2003) claim that this model represents an "empirical disruption law," 
in which different galaxies have nearly the same exponents, k ps 0.6, but timescales r* that 
differ by more than two orders of magnitude. It is therefore both mass dependent and galaxy 
dependent. Fall et al. (2009) derived analytical formulae for the mass and age distributions 
for this and other disruption models. For the LMC clusters, de Grijs & Anders (2006) 
advocate the parameter values k = 0.6 and r, = 8x 10 9 yr. As Figure 3 shows, this model 
fails by a wide margin to match the observed mass and age distributions of the LMC clusters. 
Chandar et al. (2010a, 2010c, 2011) present similar comparisons and failures for the SMC, 
M83, and M51 clusters. 

A key point illustrated by Figure 3 is that any mass dependence in the disruption rate 
will introduce features in both the mass and age distributions. These occur at M d (r) = 



10 4 (r/ r *) lAM © and r d {M) = r,(M/10 4 M Q ) fc , respectively. For M > M d (r) and r < r d (M), 



the mass function retains its initial shape and amplitude, while the age distribution remains 
flat. Conversely, for M < M d (r) and r > r d (M), the mass function becomes shallower, while 
the age distribution steepens. We note that these signatures of mass-dependent disruption 
are even more prominent in the age distribution than in the mass function and would be easy 
to detect if they were present. As we have already emphasized, there are no such features in 
the observed distributions displayed in Figures 1 and 2. 

The mass and age distributions compiled here are based on samples of clusters that are 
complete or unbiased with regard to mass and age. We caution readers that some of the 
distributions in the literature are based on samples with less suitable selection criteria. The 
completeness of the samples may be poorly known and/or vary with mass and age; clusters 
may be excluded if they are located in crowded regions, have uncertain photometry, or are 
somehow deemed to be unbound. Selection biases such as these can easily cause deviations 
from power laws in the mass and age distributions, as demonstrated explicitly by Fall et al. 
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(2009) and Chandar et al. (2010a). Stochastic fluctuations in the luminosities and colors 
of clusters are another cause of spurious features in these distributions, especially when 
relatively narrow mass and age bins are chosen (e.g., the RSG "gap" at r w (1-3) xlO 7 yr; 
Fall et al. 2005; Fouesneau et al. 2012). 

Some authors also claim to find an exponential steepening of the mass function for 
M > M c ~ 1O 5 M (Larsen 2009). These claims, however, are based on the absence of only 
a few clusters relative to extrapolated power laws and thus have little statistical significance 
(Chandar et al. 2010b). The mass functions displayed in Figure 1 do not even hint at such 
steepening. If they did, it would be appropriate to replace the power law in Equation (3) 
by a Schechter function, giving g(M,r) oc M^r 1 exp(— M/M c ). We note that this mass-age 
distribution also represents mass-independent disruption (for M c = const) because it is still 
of separable form, g(M, r) cx ?/>(M)x(r). This illustrates an important point: features in the 
mass function are not by themselves evidence for mass- dependent disruption^ 

Bastian et al. (2012) have presented new mass and age distributions for clusters in two 
fields in M83, including the one studied by Chandar et al. (2010c). We have checked that 
these distributions are mostly consistent with those derived by Chandar et al. and reproduced 
here in Figures 1 and 2. The discrepancies are relatively small and can be attributed to 
selection biases of the type mentioned above. The Bastian et al. samples are incomplete 
for young clusters of all masses (r < 10 7 yr), for older low-mass clusters (r > 10 8 yr and 
M < 1O 4 M ), and for clusters in crowded regions. When we apply these selection biases to 
the more complete Chandar et al. sample, we introduce artificial features in the mass and age 
distributions, mimicking those in the Bastian et al. distributions. We are currently analyzing 
complete samples of clusters in these and several more fields in M83 and will present the 
results in a separate paper. 



3. FORMATION AND DISRUPTION PROCESSES 

In this section, we discuss the physical processes most likely to have shaped the observed 
mass and age distributions shown in Figures 1 and 2. Before proceeding, we emphasize that 
our power-law models for i[>(M), x(t), and g(M, r) are convenient fitting formulae suggested 



2 For (3 > —2, the mass function requires an upper cutoff to prevent the total mass of young clusters in 
a galaxy from diverging. The debate is whether this cutoff occurs at a relatively low mass (M c ~ 1O 5 M ) 
and whether it has actually been detected in the data. Most of the claimed detections are based on indirect 
methods, and most ignore extinction and uncertainties in the exponents /3 and 7 of the mass and age 
distributions. 
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directly by observations and, as such, are valid irrespective of any theoretical interpretation. 
It is also worth emphasizing that we currently lack some auxiliary information needed to 
construct a definitive explanation for the mass and age distributions. There are, for example, 
no surveys of molecular clumps in galaxies beyond the Milky Way. We also lack accurate 
measurements of the internal structure of large and complete samples of extragalactic star 
clusters. Nevertheless, we have enough theoretical and observational knowledge to construct 
a partial explanation of the mass and age distributions as follows. 

The observed mass functions of molecular clouds and clumps in the Milky Way are 
power laws with exponents (3 ~ —1.7 (Shirley et al. 2003; Rosolowsky 2005; Munoz et al. 
2007; Wong et al. 2008). Because more massive clouds and clumps have longer lifetimes than 
less massive ones, they are overrepresented in the samples from which the mass functions are 
derived. Thus, the initial mass function, corrected for this effect, must be slightly steeper 
than the observed one. The result of this correction is (5 « —2.0 for clouds and clumps, nearly 
identical to the exponent for clusters (Fall et al. 2010). This is what we would expect to find 
if all the gas in the clumps (protoclusters) were converted into stars, i.e., if the efficiency of 
star formation £ were 100%. In fact, however, the efficiency is much lower: £ ~ 20%-30% 
in the dense clumps where clusters form and a factor of ~10 lower still in the larger and 
more diffuse clouds that contain the clumps (Lada & Lada 2003). Thus, to explain why the 
mass function of clusters is a power law with exponent j3 ~ —1.9, we must explain why the 
efficiency of star formation £ is independent (or nearly so) of the masses of their antecedent 
clumps. 

Stars form in a protocluster until their output of energy and momentum has expelled the 
remaining gas. This stellar feedback potentially includes radiation pressure on dust grains, 
photoionized gas pressure, protostellar outflows, main-sequence winds, and supernovae. It 
is likely that radiation pressure predominates in massive protoclusters and that supernovae, 
because they come after the other feedback processes, are relatively ineffective at removing 
gas (Krumholz & Matzner 2009; Murray et al. 2010). The condition that feedback can 
accelerate the remaining gas to the escape speed then leads to formulae for £(M) (Fall et al. 
2010). These depend on whether the feedback is closer to the energy-driven limit (with 
negligible radiative losses) or to the momentum-driven limit (with maximum losses) and 
the relation between the escape speed and the masses of the clumps or, equivalently, their 
characteristic radius-mass relation, which we approximate by a power law, R oc M a . Most 
of the feedback in the dense protoclusters is momentum driven (see Table 1 of Fall et al. 
2010), and in this limit, we have £ oc M l - 2a . 

Fall et al. (2010) showed that there is a strong correlation between the radii and masses 
of molecular clumps, with a 0.4, based on measurements of CS, C 17 0, and 1.2 mm dust 
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emission in the independent surveys by Shirley et al. (2003), Faundez et al. (2004), and 
Fontani et al. (2005). The more recent survey by Wu et al. (2010), based on measurements 
of different molecular tracers, confirms this correlation, with a ~ 0.5 (see their Figures 31- 
34). These radius-mass relations imply S(M) w constant for momentum- driven feedback. 
Thus, we have a simple physical explanation for the similarity between the mass functions 
of clumps and clusters. It remains, then, to explain the power-law form of the mass function 
of the clumps in terms of physical processes in the turbulent ISM, an important challenge 
beyond the scope of this paper. 

Following the expulsion of natal gas from a cluster, it will continue to lose mass from its 
member stars during the normal course of stellar evolution (through winds and other ejecta). 
Stellar mass loss itself removes only 10%-50% of the mass of a cluster (depending on its age, 
metallicity, and stellar IMF), but if the cluster is weakly bound and resides in a tidal field, it 
can be partially or even completely disrupted by this process (Chernoff & Weinberg 1990). 
This happens because the cluster cannot maintain a state of virial equilibrium (Fukushige 
& Heggie 1995). Clusters with low concentrations, c = log(r 4 /r c ) < 0.7 (where r t and r c are 
the tidal and core radii), are particularly vulnerable. In terms of the dimensionless central 
potential, this condition is equivalent to W < 0.3. We expect such weakly bound clusters 
to be numerous as a result of the earlier expulsion of gas by stellar feedback. 

The effect of this disruption process on the mass and age distributions of the clusters 
depends primarily on their concentrations — for clusters with the same stellar IMF, until 
two-body evaporation becomes important (Chernoff & Weinberg 1990; Fukushige & Heggie 
1995). In particular, the shape of the mass function will be preserved, while its amplitude 
declines, so long as c is not correlated with M. This is a reasonable expectation because 
the star formation efficiency £, which determines the relative weakening of the gravitational 
potential of the clusters by stellar feedback, is independent of M. The top panel of Figure 4 
shows c plotted against M for LMC and SMC clusters from fits of King profiles to Hubble 
Space Telescope (HST) images by McLaughlin & van der Marel (2005). Evidently, there is 
no correlation between c and M, consistent with the fact that the mass function has the 
same power-law shape at all ages. This result is suggestive rather than definitive, however, 
because it involves only a small fraction of clusters in the LMC and SMC@ 

Star clusters can also be disrupted by tidal interactions with passing molecular clouds. 
Spitzer (1958) proposed this mechanism to account for the scarcity of clusters older than 



3 The McLaughlin & van der Marel (2005) sample includes most of the LMC and SMC clusters with 
structural parameters derived from HST images. We have compared the values of r c , r^, and r t derived from 
HST and ground-based images of the same clusters and find large discrepancies, indicating that the latter 
are not accurate enough for our purposes. 
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~10 9 yr in the galactic disk. We follow the comprehensive treatment by Binney & Tremaine 
(2008) and distinguish two regimes: catastrophic, in which clusters are disrupted suddenly 
by a single strong encounter, and diffusive, in which clusters are disrupted gradually by a 
series of weak encounters. We denote the mass and half-mass radius of the clusters by M 
and Th and the corresponding quantities for the perturbers (molecular clouds) by M p and 
r\ p - The characteristic internal density of the clusters is then given by ph = 3M/(87rr^). 
Furthermore, we denote the mean number of perturbers per unit volume by n p and the RMS 
dispersion of relative velocities between the clusters and perturbers by a v . 

The timescale for disruption of clusters by tidal interactions depends on these parameters 
in the two regimes as follows: 

1/2 



Td ^ la — (catastrophic regime), (6) 

2 

Td ^ ^72^~ (diffusive regime). (7) 

In both regimes, Td depends on the properties of the clusters only through p h . Thus, this 
mechanism will not change the shape of the mass function of the clusters unless there is a 
correlation between ph and M. 

The middle and bottom panels of Figure 4 show the half-mass radius and density plotted 
against mass for LMC and SMC clusters with HST images analyzed by McLaughlin & 
van der Marel (2005). Evidently, there is a positive correlation between r h and M but 
none between ph and M, consistent again with the mass function preserving its power-law 
shape. As before, this result is suggestive rather than definitive, because it is based on an 
incomplete sample of LMC and SMC clusters. A correlation between p^ and M would be 
surprising, however, because the tidal field of the host galaxy tends to impose the same mean 
density, p t = 3M/(47rrj i ) = 2(rf l /r t ) 3 pf l , on all clusters at the same galactocentric distance, 
independent of M. This constraint permits a correlation between ph and M only if there is 
a compensating correlation between rh/r t ~ 0.7(r c /r t ) 1 / 2 and M and hence between c and 
M. As noted above, such a correlation is neither expected nor observed. 

The disruption timescale in the catastrophic regime depends on the properties of the 
molecular clouds (perturbers) only through their mean smoothed-out density, p p = M p n p . 
Binney & Tremaine (2008) estimate that open clusters in the solar neighborhood marginally 
satisfy the criterion for disruption in the catastrophic regime and have ~ 3 x 10 8 yr. This 
local estimate of r d should also apply to other galaxies, after scaling inversely by p p , if the 
internal densities ph of the clusters are known or assumed to be the same as those in the 
solar neighborhood. The situation is more complicated in the diffusive regime because Td 
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then also depends on M p , r^ p , and a V: quantities that are much harder to determine from 
observations than p p . In this case, Equation (7), which was derived for identical perturbers, 
should be revised so that the disruption rate 1/r^ is a weighted sum over terms with different 
M p and r^, based on their frequency of occurrence in the population of molecular clouds. 

Over longer times, low-mass clusters are disrupted by stellar escape driven by inter- 
nal two-body relaxation ("evaporation"). This process, unlike the others discussed here, 
inevitably causes a bend in the mass function. The mass at which the bend occurs is age 
dependent, reaching M ~ 1O 5 M at r ~ 10 10 yr, which provides a natural explanation for 
the observed turnover in the mass function of old globular clusters (Fall & Zhang 2001; 
McLaughlin & Fall 2008 and references therein). However, we do not expect to observe this 
feature in the mass functions displayed in Figure 1, because in this case, the clusters are too 
massive and too young to have experienced much two-body evaporation (as demonstrated 
quantitatively in the papers from which the mass functions were taken). 

There are some reasons for expecting the mass and age distributions of clusters to be 
similar in different galaxies. First, the mass functions of molecular clouds are observed to be 
broadly similar (Rosolowsky 2005; Blitz et al. 2007; Fukui et al. 2008; Wong et al. 2011)0 
If the mass functions of the clumps, which have yet to be surveyed outside the Milky Way, 
are also similar, then clusters in different galaxies would form with similar initial conditions. 
Second, the disruption by feedback and subsequent stellar mass loss, processes internal to 
the clusters, is expected to operate similarly in different galaxies. The rate of disruption by 
passing molecular clouds may differ among galaxies, depending on the mean smoothed-out 
density p p and other factors (see Equations (6) and (7)), but this effect may be too weak to be 
detectable over much of the accessible range of ages. For example, the estimated disruption 
timescale for clusters in the solar neighborhood, ~ 3 x 10 8 yr (Binney & Tremaine 2008), 
exceeds the ages of all the clusters used to construct the age distribution plotted in Figure 2. 

Some of the disruption processes discussed here are likely influenced by the tidal field 
of the host galaxy. Since this has a strong inverse dependence on distance from the galactic 
center, one might expect the mass and age distributions of the clusters to exhibit some spatial 
variation. Unfortunately, it is difficult to make reliable predictions of this effect because it 
depends on the initial concentrations of the clusters, which in turn are governed by the 
expulsion of gas by stellar feedback. In particular, it is not clear whether there would be 
much of a correlation between the initial concentrations and the locations of clusters. In any 



4 Some differences have been reported in the mass functions of molecular clouds in different galaxies. 
These comparisons are based on data from different instruments, analyzed in different ways by different 
authors, and pertain to different density thresholds and mass ranges. Thus, the reported differences probably 
exaggerate the real differences to some extent. 
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case, it is important to note that the characteristic tidal fields within galaxies — measured, 
for example, at their effective radii — have relatively weak variations among galaxies (Fall & 
Zhang 2001; McLaughlin & Fall 2008). Thus, although the disruption of clusters may be 
influenced by local tidal fields, we do not expect this to cause major differences in the mass 
and age distributions when averaged over whole galaxies. 

4. DISCUSSION 

We have reexamined the mass and age distributions of star clusters in the Milky Way, 
LMC, SMC, M83, M51, and Antennae galaxies. These distributions are well represented 
by power laws: ip(M) oc dN/dM oc M 13 with j3 ps —1.9 and x( T ) °c dN/dr oc r 7 with 
7 w —0.8. Furthermore, the mass and age distributions are approximately independent 
of each other: g(M,r) oc ip(M)x(j). Since the rates of star and hence cluster formation 
have varied relatively slowly, the more rapid decline of the age distributions must be the 
result of cluster disruption. We find no evidence for mass-dependent disruption; indeed, 
simple models with oc M os fail by wide margins to match the observed mass and age 
distributions. 

The mass functions of the clusters are similar to those of the molecular clouds and clumps 
in which they form, requiring a nearly constant (but small) efficiency of star formation. This 
is an expected consequence of stellar feedback, given the observed radius-mass relation of the 
clumps (R oc M a with a ~ 0.5). The clusters are also disrupted by subsequent stellar mass 
loss and passing molecular clouds. These processes preserve the power-law shape of the mass 
function unless there are correlations between the concentration c and M or between the 
internal density ph and M (until two-body evaporation becomes important). The available 
data, while limited, reveal no such correlations. 

The main message here is the similarity among the mass functions and age distributions 
of clusters in different galaxies, strikingly evident in Figures 1 and 2. The literature on 
this subject includes some contradictory claims: for prominent features in the mass and 
age distributions, for strong dependencies on mass or differences among galaxies. We have 
shown that selection biases are responsible for several of these claims and we suspect they are 
responsible for the others. It is hard to see how selection biases and/or analysis errors would 
make the mass and age distributions appear more similar than they really are, whereas it is 
easy to see how such problems could introduce spurious differences. In any case, since the 
distributions presented here are based on complete or unbiased samples, we are confident 
that the similarity they exhibit is real. 
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Nevertheless, we do not expect the mass and age distributions to be exactly the same 
in all parts of all galaxies. As we have emphasized here and in our previous papers, we 
expect minor differences in the age distributions as a consequence of different histories of 
cluster formation and/or disruption by passing molecular clouds. We would also not be 
surprised to find minor differences in the mass functions of molecular clumps and hence 
those of young clusters. Indeed, we find minor differences in the observed distributions; the 
best-fit exponents have dispersions ap = 0.15 and <r 7 = 0.18. The reality of these differences 
remains questionable, however, because o$ and <r 7 are close to the true uncertainties in the 
exponents, ~ e 7 ~ 0.1-0.2. 

We also expect any differences in the mass and age distributions to depend on spatial 
scale. In small regions (< few x 10 2 pc), there may be large variations in the formation 
and disruption rates, and hence in the mass and age distributions. As more of these small 
regions are combined into larger ones, the variations will average out, and the differences in 
the mass and age distributions will diminish. Figures 1 and 2 show the important result that 
these differences are negligible or barely detectable on the scale of whole galaxies. Evidently, 
the similarity among the mass and age distributions overwhelms such differences. This is 
the sense in which we consider the distributions to be "universal" or "quasi-universal." 

Finally, we mention several future studies that would help to advance this subject. On 
the observational side, it would be interesting to survey the molecular clumps in nearby 
galaxies, to derive their mass functions, radius-mass relations, and other properties, for 
comparisons with those in the Milky Way. Another important goal is to measure the internal 
structure of clusters from EST images (especially c and ph) in large unbiased samples in 
several nearby galaxies. From our experience, Ha measurements are necessary for accurate 
determinations of the mass and age distributions of young clusters (r < 3 x 10 7 yr) and 
should be included in all future studies. On the theoretical side, it would be interesting to 
explore further the effects of stellar mass loss on the evolution, stability, and disruption of 
weakly-bound clusters in tidal fields. 
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Fig. 1. — Mass functions of star clusters in different age intervals in different galaxies (as 
indicated). These have been adapted from the references given in the text. The absolute 
normalizations of the mass functions are arbitrary, but the relative normalizations within 
each panel are preserved. The lines show power laws, dN/dM oc M 13 , with the best-fit 
exponents listed in Table 1. Note that these are all close to f3 — —1.9. 
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Fig. 2. — Age distributions of star clusters in different mass intervals in different galaxies 
(as indicated). These have been adapted from the references given in the text. The absolute 
normalizations of the age distributions are arbitrary, but the relative normalizations within 
each panel are preserved. The vertical spacing between the age distributions depends on 
the adopted mass intervals, which differ among the galaxies for practical reasons (distance, 
limiting magnitude, sample size). The lines show power laws, dN/dr oc r 7 , with the best-fit 
exponents listed in Table 2. Note that these are all close to 7 = —0.8. 
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Fig. 3. — Mass functions and age distributions of star clusters in the LMC (left-hand and 
right-hand panels, respectively). The data points are the same as in Figures 1 and 2. The 
curves are based on a simple model with gradual mass-dependent disruption, computed from 
Equation (12) of Fall et al. (2009) with f3 — -2, k — 0.6, and either r* = 8 Gyr (top panels) 
or = 1 Gyr (bottom panels). The three curves in each panel correspond to the same age 
and mass intervals as the data points. Note that the models are poor fits to all the data, 
except the mass function at the youngest ages (r < 10 7 yr). 
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Fig. 4. — Concentration, half-mass radius, and half-mass density plotted against mass for 
LMC clusters (circles) and SMC clusters (triangles) in two intervals of age (as indicated). 
All quantities were determined by McLaughlin & van der Marel (2005): c, r h) and ph from 
fits of King models to HST images and M from total luminosities and mass-to-light ratios 
inferred from integrated colors and stellar population models. Note that there is a positive 
correlation between r\ and M but none between c and M or between ph and M. 
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Table 1. Exponents of Mass Functions 



Galaxy 


Age Interval 
log(r/yr) 


P 1 


Milky Way 


<6.5 


-1.82 ±0.01 


LMC 


6-7 


-1.70 ±0.02 


LMC 


7-8 


-2.02 ±0.23 


LMC 


8-9 


-1.87 ±0.07 


SMC 


6-7 


-1.87 ±0.25 


SMC 


8-9 


-1.84 ±0.14 


M83 


6-7 


-1.91 ±0.10 


M83 


7-8 


-1.86 ±0.13 


M83 


8-8.6 


-2.06 ±0.13 


M51 


8-8.6 


-2.24 ±0.08 


Antennae 


6-7 


-2.14 ±0.03 


Antennae 


7-8 


-2.03 ±0.07 



^east-squares fits to \og(dN/dM) = 
(5 logM ± const. 
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Table 2. Exponents of Age Distributions 



Galaxy 


Mass Interval 
log(M/M ) 


7 1 


Milky Way 


> 2.2 


-0.97 ±0.16 


LMC 


2.8-3.4 


-0.54 ±0.10 


LMC 


3.4-4.0 


-0.68 ±0.07 


LMC 


4.0-4.8 


-0.75 ±0.07 


SMC 


3.0-3.5 


-0.73 ±0.17 


SMC 


3.5-4.5 


-0.70 ±0.16 


M83 


3.5-4.0 


-0.94 ±0.09 


M83 


> 4.0 


-0.95 ±0.09 


Antennae 


4.3-5.0 


-1.04 ±0.30 


Antennae 


> 5.0 


-1.05 ±0.11 



^east-squares fits to \og(dN/dr) = 7 log r 
+ const. 



